Spin polarization and magneto-luminescence of confined electron-hole systems 
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A BCS-like variational wave-function, which is exact in the infinite field limit, is used to study 
the interplay among Zeeman energies, lateral confinement and particle correlations induced by the 
Coulomb interactions in strongly pumped neutral quantum dots. Band mixing effects are partially 
incorporated by means of field-dependent masses and g factors. The spin polarization and the 
magneto-luminescence are computed as functions of the number of electron-hole pairs present in the 
dot and the applied magnetic field. 
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I. INTRODUCTION 

In the present paper, we study the spin polarization 
and the magneto-luminescence of a neutral, medium-size 
quantum dot (qdot) subjected to a strong (pulsed) laser 
pumping and a strong magnetic field. There are many 
good reasons to study the properties of this system. 

On one hand, recent magneto-tunneling experimentd 
have stated very clearly the role of Zeeman energies, lat- 
eral confinement and Coulomb repulsion in the spin po- 
larization of a qdot filled with N electrons. At "filling 
factors" 2 > ly ^ 1, ground-state rearrangements lead to 
significant oscillations of the conductance peak positions 
as a function of the magnetic field. The situation could 
be even more interesting for electron-hole (e-h) clusters, 
where additional e-h correlations come into play. Traces 
of spin rearrangements shall be seen in the luminescence 
and other optical properties. 

On the other hand, extensive measurements of quan- 
tum well (qwell) luminescence exist for different e-h. den- 
sities (laser excitation power) and Bolarizationscl. A 
BCS-based theory has been proposedB. Very high mag- 
netic fields, up to 60 T, have been applied mainly to 
low-density systems, and the results have been inter- 
preted iB-terms of isolated neutral and negatively charged 
excitongj. Analogous experiments in qdots are lack- 
ing, but the available experimental resources are enough 
to create a high population of excitons in a medium- 
size qdot in a strong magnetic field. 6-exciton lumines- 
cence has been undoubtedly observed in single small InAs 
quantum dots at _B = CO. In qwells, very high e-h den- 
sities {iQi"^ pairs/cm^) have been achieved with pulsed 
pumpingQ. 

Our model quantum dot is made up from a symmet- 
rical GaAs-AlGaAs quantum well, 8.5 nm wide in the 
growth direction. Stress is supposed to induce a lateral 
confinement, which is described by a parabolic potential 
with hH. « 1.2 meV. An electron- hole system is created 
in the dot by a strong 5 ps pulsed laser, as it is experi- 
mentally done for example in Ref|7[ The mean lifetime of 
this system is around hundred of ps or even longer, i. e. 
a time scale much greater than the characteristic times 
(« 1 ps) to reach equilibriumB. Thus, at very low tem- 



peratures we end up with an "stable" A^-pair e-h cluster 
in its ground state. 

Results for spin- up and -down densities, hole and elec- 
tron spin polarizations and for the position and magni- 
tude of the coherent luminescence peak as functions of 
the magnetic field and of the mean number of excitons 
are presented below. The theoretical framework used 
is a BCS-like wave function, corrected against particle 
number navi-ponservation by means of a Lipkin-Nogami 
procedurcBllj. This wave function is able to reproduce 
the e}i ] pO |Cted "Bose condensed" state in the B ^ oo 
limitllJ'tS'EJ. A big basis of one-particle states is used, 
which includes up to 3 Landau levels and 202 states per 
Landau level. We consider systems with up to 40 e-h 
pairs. 

The plan of the paper is as follows. In Section II, the 
model to be employed is described in details. The basics 
of our theoretical scheme are summarized in Section III. 
In Section IV, the main results are presented and dis- 
cussed. Concluding remarks are given in the last section. 



II. THE MODEL 

We consider a system of N electrons and N holes con- 
fined in a quasi two-dimensional quantum dot, and in the 
presence of a perpendicular magnetic field. As mentioned 
above, the qdot is made up from a 85 A-wide symmet- 
ric qwell. A parabolic potential confines the motion of 
the particles in the plane perpendicular to the grow di- 
rection. The first qwell sub-band approximation is used, 
i. e. the confining energies along the z direction are 



written as 



— P 



ith fc. 



(h) 



= 1. Notice 



^e{h) 2m 

that for the 85 A-width GaAs well, the gap to the sec- 
ond qwell state is Ai;| « 210 meV, and A£;^ « 39 meV 
for electrons and holes respectively, whereas the typical 
Coulomb energy is Ecouiomb ~ 'i.l%^J B\T] meV. This 
model is a common theoretical framework the study 
of strained or self-assembled quantum dotsliJ. By using 

the symmetric gauge, A[r) — dimensionless coordi- 
nates, r — > Ibt, with Ib = ^/ificj eU, and using Landau 
level (LL) states = jrii, m.^, s^)} (jii = 0, 1, • • ■ ; rrii = 



2 



— cxD, • • ■ , —1, 0, 1, • ■ • , oo; and Si = ±1/2 are the radial, 
angular momentum and spin quantum numbers, respec- 
tively) as set of one-particle states, we can write the 
Hamiltonian in second quantization as: 



hLuf. 



Zeeman [ \ 
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(1) 



where rrieirnh) are the electron and hole effective masses, 
ri is the dot confining frequency, and e is the dielectric 
constant. Egap is the gap separation between the con- 



duction and valence bands, -E^^^™"" 



e(h) 



ton, s 



with 5e(/i) the effective g-factors, /i is the Bohr magne- 
''^^ are the z-components of the ith particle spin, 
2ni + jm^l ± -|- 1 are the LL energies for elec- 
trons (holes), — eB/me(^ii)C is the electron (hole) 
cyclotron frequency, and ei{hi), e-(/i-) are the electron 
(hole) destruction and creation operators. Convention- 
ally, we write Sh = ±1/2 for the two branches of the 
heavy hole sub-band. To the electronic — —3/2, for 
example, we ascribe Sh = 1/2. 

The effective parameters entering the Hamiltonian 
(masses and g factors) are magnetic field- and width- 
dependent magnitudes to approximately account for 
band mixing. For the 85 A welL we fitted the experi- 
mental in-plane heavy hole massc3, thus obtaining: 



,85 



A 



iB) = 



0.17, 

0.17+0.0168B 
1+0.023B ' 



for B < 10 T 
for S > 10 T. 



(2) 



Experimentally, the dependence of the electron g- 
factor on , well ..width and magnetic field is well 
determinedE^llZHiSo'EJ, and in our case we have: 



|With the field, as in the electron case, and fitted it to 
The result is 



.85. 



A 



(B) = -0.05 B[T]. 



(4) 



Notice that in the i? — » limit, gt vanish because the 
exciton (X) g-factor, defined as gx = ge + gh, ^sd the 
electron g-factor are equal in the 85 A width welllla. 

With this parameterization, ge changes sign at Be ~ 32 
T. Consequences of this facts are discussed below. 



III. THE BCS SCHEME 

The BCS approach has been successfully applied by 
Paquet et. -al. in the study of the two-dimensional (2D) 
e-h systemtj, and by Feraandez-Rossier and Tejedor to 
the exciton gas in a qwellQ. We used it previousljpp the 
study of finite e-h systems at zero magnetic fieldEHl. We 
employ the Lipkin-Nogami (LN) scheme to avoid particle 
number non-conservation in a finite system. The BCS 
wave function is given by 



\BCS)n = n + ^^hl) \0)h\0). 



(5) 



where \0)h{\0)e) are the respective electron (hole) vac- 
uum. The subscript "iV" in the BCS function means 
that the average number of pairs is just N, i.e. (TV) = 
{BCS\N\BCS)n = N. Ui and Vi are the variational 
parameters. They fulfill the normalization conditions 
^i+'i'i= li The hole state \t) which is paired with 
the electron state \i) = \ni,mi,Si) is different for the 
two magnetic field regimes, reflecting the ground state 
spin structure, as it is show in Fig. |^. For B < 32 T, 
electrons and holes with opposite spins are paired, then 
|z) = jrij, — mj, —Si). Whereas for i? > 32 T, the energy 
is minimized when the \i) electronic state is paired with 
\i) = |?ij,-mi,Si). 

The total angular momentum (projection onto the z 
axis), corresponding to \BCS)n is zero because the an- 
gular momentum of each pair is zero. The total electron 
or hole spin, however, depend on the populations of spin- 
up and down components. 

The detailed—description of the LN method can be 
found in Refs.Q'EB. For completeness, we sketch the main 
results. The LN estimate for the ground state energy is 
given by: 



ge{B) = -0.1667 + 0.0052 B[T]. (3) 

The dependence of the hole g-factor on the mag- 
netic field for iigh B values, however, is not prop- 
erly determinedESIij. Here we assume a linear behavior 



Eln = Ebcs - 2Xi{{N) -N)- X2{{N^) - N^), (6) 

where Ebcs is the expectation value of 

the effective Hamiltonian H — iJ — 
N{Egap - El{ke = 1) - Elikh = 1)} in \BCS)n- 



3 



B < B„ 



B> B 



Conduction 
Band 



Valence 
Band 




5" 
E 



30 35 
Magnetic Field [T] 




FIG. 1: Schematics of the spin sub-band structure for the two 

field regimes. The dotted lines show the paired electron and 
hole states in the BCS wave function. The interband optical 
transitions for the two circular polarizations of light cr* are 
indicated. 
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Notice that in the second term, the sum runs over hole 
states \i). ef = £j = 2ni + Imjl + + 1 are the one- 



particle energies, and e\ 



OA 



.OA 



2ni + m, + 1 



{i\r \i) is the expectation value of the harmonic poten- 
tial. The mean value of the number of pairs is 

i i 

= E^l (8) 

i 

The extrema conditions can be written in the standard 
form of gap equations 



(9) 



FIG. 2: Dependence of Eopt — E^ark with the magnetic field. 
In figures 2-5, results for the iV = 40 exciton system are 
presented. 



where the Hartree-Fock energies are given by 



HF 



-{hujl + J^){2n, + \mi\ + mi + 1) 

1 fhnl hnl \ , , 
o -f + -7r (2«^ + KI + i) 



UJ" 



u 1 
+ -^Big^Si - ghS^) - ^("1 j^l") 



els 



E {^J\^\J^W,-MN-v}), (10) 



and we have used the usual BCS parameterization 



(11) 



The chemical potential /x = Ai -|- A2/2 was introduced to 
fix the particle number, and A2 is determined in the LN 
scheme as: 



A2 = {{HN^)[n^-{N^))+{HN)[{N^)-{NN'^)) 
+ {H) (^{Ny - {NN^)^ } {(iV2)3 + {N^f 

+ N^N^) - {N^) (2N{N^) + {N^)^Y^ , (12) 



where the expectation values (• • •) arc taken in the 
\BCS)n state. The resulting equations are solved it- 
eratively up to a precision of 10~^^ in ef^. Calculations 
were performed for < N < AQ pairs and 606 one- 
particle LL states. 
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IV. POLARIZATION AND 
MAGNETO-LUMINESCENCE 

As shown below, the spin polarization of the electronic 
subsystem increases as the magnetic field is increased. At 
B ~ 32 T, the electronic g-factor changes sign according 
to Eq. 1^, leading to a rearrangement of electron-hole 
pairing in the absolute ground state. The kind of pairing 
minimizing the energy is represented in Fig. |l| by dashed 
lines. For B > 32 T, the ground state is dark, and it is not 
even clear whether it can be reached by means of light 
excitation followed by spin relaxation processes. Thus, 
besides the ground state, for i? > 32 T, we compute also 
the lowest BCS state with (t+ and a~ excitons. Below, 
we present results for the TV = 40 system, obtained with 
three LLs and 202 states per level, i. e. a total of 606 
one-particle states. 

As mentioned above, two BCS functions may be con- 
structed. One in which optical excitons are formed, and 
a second one in which dark excitons are formed. The dif- 
ference Eopt — Edark is drawn in Fig. |^ as a function of 
B, showing that the dark state becomes the ground state 
when the value i? « 32 T is crossed and the electronic 
sub-bands are re-ordered. The absence of efficient spin 
relaxation mechanisms may, however, prevent the actual 
ground state to be occupied. 

It is interesting to note that the difference Eopt— Edark, 
in the magnetic field interval shown in Fig. ||is very close 
to the difference between the electronic Zeeman energies. 
A simple qualitative picture can be offered for the un- 
derstanding of this and the next figures. The properties 
of the system are roughly determined by the holes be- 
cause of the competition among the hole Zeeman energies 
and the total harmonic and Coulomb energies. The hole 
occupations are thus very similar for optical and dark 
states. The form of e-h pairing provides the "fine struc- 
ture" . That is, the minimization of the energy leads to a 
definite pairing. 

Scaled spin-up and -down densities, obtained from 



N 



J2 (t>'i{f^^',if^N{BCS\ele,\BCS) 

(13) 



i={ni,mi,{T ,1)} 



are shown in Fig. |. At B = 40 T, we show the 
ground (dark) and optically-active excited-state densi- 
ties, which are almost inverted in agreement with the ar- 
gument given above. Notice that there are only four ex- 
cess spin-up electrons at S = 20 T. These small net polar- 
izations for high magnetic fields are related to the attrac- 
tive character of the e-h interaction. Unlike pure electron 
systems, small-radius orbits maximize e-h attraction, and 
the competition between Zeeman and Coulomb energies 
starts at higher fields. 
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FIG. 3: Spin-up and down densities, (a) Ground state at 
B = 20 T, (b) Ground (dark) state at B = 40 T, and (c) 
Excited (optically active) state at B = 40 T. 



The "hole dominance picture" leads to changes in the 
polarization, as B is increased, through the reconstruc- 
tion of the droplet edge- in a way very similar to electrons 
near filling factor onec2l. This fact is illustrated in Fig. 
^, where the difference between ground-state spin-up and 
-down densities for different values of B are shown. It is 
evident that polarized densities differ mainly at the edge, 
and that changes in the polarization are more significant 
at the edge. 

The following two figures, Figs. ^ and ^, contain the 
main results of our paper. Electron and hole spin polar- 
izations are drawn in Fig. ^ as a function of B. Solid 
lines refer to the ground state, while the dashed lines for 
_B > 32 T refer to the optically active state. Notice that 
even at a high field value like 45 T, the electronic polar- 
ization is only 70 %. Notice also the change in sign of the 
ground-state electronic spin a,i B = Be- The inset shows 
the total ground-state spin squared, computed from 
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FIG. 4: Difference between ground-state spin-up and -down 
densities for various magnetic field values. 



{S!)BCS = :rN- 



E2 2 



(14) 

The total (coherent) magnetolumincscence intensity 
for both cr"*" and a~ polarizations is presented in Fig. 
We compute it for <7~ polarization, for example, from the 
expression 
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FIG. 5: (a) Electron and hole spin polarizations. For B > 32 
T, both ground-state (solid line) and excited-state (dashed 
line) properties are drawn. The inset shows the total electron 
spin squared, (b) Luminescence intensities for both a~ and 
cr"*" polarizations. 
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where |/) is a basis of — 1 particle states , a = (n, m) 
is a composed index and, P„- — Ea^aT^ai i'^" 
terband dipole transition operator for the a~ circularly 
polarized light. Notice that Itotai is the integrated lumi- 
nescence, corresponding to the transition from the given 
initial BCS state to any final state. The convention for 
solid and dashed lines is the same as in Fig. ^ja. It shall 
be stressed that the degree of polarization, defined from 



TO _ ja^ 

total total 

ja- I rcr+ ' 
^total ^ ^total 



(16) 



follows very well the behaviour of (iSz), i. e. the difference 
between the occupation of spin- up and -down sub-bands. 
This polarization is nearly 10 % at 20 T, and around 70 
% at 45 T. 

Finally, in Fig. ^ we show results for the position of 
the luminescence line and the intensities as functions of 



the numbers of pairs in the dot. In our computations, the 
energy is corrected against non-conservation of the total 
number of excitons, N . Conservation of the number of 
(T+ and a~ excitons is not properly taken into account. 
Thus, we can not exactly compute the position of the 
(7+ and (7^ lines. In place of it, we show in Fig. ^ the 
difference /\En = Eln{N) - Eln{N - 1) relative to the 
exciton a~ line. Ex- The following interesting properties 
can be noticed in thia.figure: a) A blueshift as the number 
of excitons is raisedQ. It is around 0.2 meV/exciton at 
N = 20, and 0.15 meV/exciton at TV = 30, and b) An 
apparent minimum of each curve at B around 30 T, i. e. 
at Be- On the other hand, the intensity (Fig. ^d) shows 
an increase with N for both polarizations, as one would 
expect from coherent emission. 



V. CONCLUDING REMARKS 

We have computed the spin polarization and the lu- 
minescence of a quantum dot in which a mean number 
of electron-hole pairs, N , have been created by a laser 
pulse. Band mixing effects were approximately taken 
into account by means of well width- and magnetic field- 
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FIG. 6: Dependence of the position of the luminescence line 
(see main text) and the intensities on the number of pairs in 
the dot. Ex is the a~ exciton energy. 



dependent masses and g-factors. For the model under 
study, the electron g-factor vanishes at Be « 32 T. It 
means that, for magnetic field values around the elec- 
tron polarization, and thus the ratio of intensities given 
by formula ( |l6| ) , is determined as a result of the interplay 
among Coulomb, confinement and hole Zeeman energies. 
The net polarization is only 70 % at i? = 45 T because 
of the attractive electron-hole interaction. 

The general features found in our calculations, i. e. 
relatively small polarizations even at high magnetic field 
values, blueshift of the luminescence lines with an in- 
crease of the laser power, etc seem to be not related to 
the specific parametrization used for carrier masses and 
(/-factors. 

The developed computational scheme may be applied 
to many other interesting situations, from which two of 
them may be distinguished. The first is the stationary 
regime, in which constant populations of and ex- 
citons arise as a result of appropriate pumping, recombi- 
nation and spin-flip processesd. The second is the study 
of the effects of hyperfine interactions between nuclear 
and electronic spins on the positipa of the recombination 
lines, known as Overhauser shift£j. Research along both 
directions is in progress. 
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